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There is a certain consensus that the very fast growth of the relaxation time r occurring in glass- 
forming liquids on lowering the temperature must be due to the thermally activated rearrangement 
of correlated regions of growing size. Even though measuring the size of these regions has defied 
scientists for a while, there is indeed recent evidence of a growing correlation length ^ in glass- formers. 
If we use Arrhenius law and make the mild assumption that the free-energy barrier to rearrangement 
scales as some power t/j of the size of the correlated regions, we obtain a relationship between time 
and length, T log t ^ . According to both the Adam-Gibbs and the Random First Order theory 
the correlation length grows as ^ ~ (T — Tt)~^^^''~^\ even though the two theories disagree on 
the value of 6. Therefore, the super- Arrhenius growth of the relaxation time with the temperature 
is regulated by the two exponents ip and through the relationship Tlogr ~ (T — Jj.)~*/(''~*'. 
Despite a few theoretical speculations, up to now there has been no experimental determination of 
these two exponents. Here we measure them numerically in a model glass-former, finding = 1 
and 6 — 2. Surprisingly, even though the values we found disagree with most previous theoretical 
suggestions, they give back the well-known VFT law for the relaxation time, Tlogr ~ (T ~ Tk)~^ . 



I. INTRODUCTION 

One of the most interesting open problems in the 
physics of glass-forming liquids is how the slowing down 
of the dynamics is related to the presence of a growing 
correlation length. Even though the very steep growth of 
the relaxation time on lowering the temperature is per- 
haps the experimentally most conspicuous trait of glassy 
systems, the existence of an associated growing length- 
scale has been a matter of pure speculation for quite 
a long time. Yet, the qualitative idea that time grows 
because relaxation must proceed through the rearrange- 
ment of larger and larger correlated regions, is a sound 
one, so that many analytical, numerical and experimen- 
tal efforts have been devoted in the last fifteen years to 
detect correlated regions of growing size in supercooled 
liquids. 

The first breakthrough was to discover the existence 
of dynamical correlation length ^d- This is the spatial 
span of the correlation of the mobility of particles and it 
is directly related to the collective dynamical rearrange- 
ments of the system More recently, a static correla- 
tion length ^ has been measured, by studying how deeply 
amorphous boundary conditions penetrate within a sys- 
tem @, [2II. Both lengthscales grow when decreasing T, 
and even though the relationship between dynamic and 
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static correlation length is still under investigation (sj, 
their existence strengthens the idea of a link between 
time and length in glass-formers. 

In principle, to uncover the formal nature of such link 
one simply needs to compare the dependence of both re- 
laxation time t{T) and correlation length £,{T) on the 
temperature, and infer the law t = ''"(1^). In practice, 
data on the correlation length are still too scarce to pur- 
sue this indirect (parametric) road reliably. Hence we 
must use a direct method to link length and time. 

The super- Arrhenius growth of the relaxation time in 
fragile systems suggests that the free-energy barrier A to 
relaxation must grow as well on lowering T. It is therefore 
natural to assume that this barrier grows because the size 
of the regions to be rearranged, i.e. the correlation length 
^, does. A mild assumption is that the following, 

A^e, (1) 

so that the Arrhenius law of activation gives 

r~exp(^^y (2) 

Of course, this crucial equation is of little help as long 
the value of the exponent ip remains unknown, but this 
is, unfortunately, the current state of affairs. 

To work out a more explicit relationship between time 
and temperature, we need to know how precisely the 
correlation length increases with decreasing T. This is, 
however, one of the most disputed and open problems 
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currently in the physics of glassy systems, so that try- 
ing to follow a common path from now on is hopeless. 
We choose to move in the framework of two theoretical 
schemes that, despite many conceptual differences, share 
some common ground, namely the Adam-Gibbs (AG) 
and the Random First Order (RFOT) theories. The 
two common ideas of AG and RFOT are: 1) there exists 
a static correlation length, whose growth is responsible 
for the growth of the relaxation time, and 2) the static 
correlation length grows because the configurational en- 
tropy Sc decreases. Although the two theories are really 
quite different in explaining how and why the correlation 
length is connected to the configurational entropy, they 
both predict a power-law link. 



(3) 



where d is the space dimension. The exponent 9 is zero 
according to AG (in fact, it is not even introduced in the 
theory). 



0AG = 0, 



(4) 



whereas it has a crucial role within RFOT, where it regu- 
lates how the interfacial energy of amorphous excitations 
grows with the size of the excitations. The exponent 9 is 
subject to the constraint 9 < d — 1. The value 



7RFOT 



d/2 



(5) 



for 9 was proposed in from renormalization group ar- 
guments. Independently of the exact value of 9, we see 
that as long as 9 > 0, the growth of ^ with decreasing 
configurational entropy predicted by RFOT is sharper 
than in AG. 

From the behaviour of the (extrapolated) excess en- 
tropy of the supercooled liquid with respect to the that 
of the crystal, it is possible to see that the configurational 
entropy scales almost linearly with T, 
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where Tk is the so-called Kauzmann's temperature, or 
entropy crisis point. In this way, the relationship between 
length and temperatures in the context of the AG and 
RFOT schemes becomes 
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T-Tk 

and that between time and temperature, 

A 

exp 



(7) 



(8) 



where A is a factor weakly dependent on temperature. 

Hence, we see that the two exponents ip and 9 con- 
tain all the relevant information about the relationship 



between time, length and temperature in supercooled liq- 
uids. However, little is known about them. In this work 
we make a numerical determination of both exponents; 
this is what we mean by a direct method to link correla- 
tion length and relaxation time. 

According to AG 5] , the barrier scales like the number 
of particles involved in the rearrangement, so that 



V'AG 



(9) 



whereas in RFOT [J], the exponent tp is fixed by a nu- 
cleation mechanism to be equal to the interfacial energy 
exponent. 



V'RFOT = 9. 



(10) 



Notably, despite predicting quite different values of ip 
and 9, in both AG and RFOT relation ^ reduces to the 
classic Vogel-Fulcher-Tamman law of relaxation, 



TvFT ^ exp 



A 



T-Tk 



(11) 



Here we find for ip and 9 values different both from AG 
and RFOT, but surprisingly still consistent with a plain 
VFT relaxation form. 



II. NUMERICAL PROTOCOL FOR THE STUDY 
OF AMORPHOUS EXCITATIONS 

To estimate the two exponents 9 and ^ we must mea- 
sure the energy of the amorphous excitations that sponta- 
neously form and relax at equilibrium. In particular, the 
exponent 9 regulates the growth of the interface cost with 
the size of the excitation, while the exponent ip rules the 
dependence of the relaxation time on the size of cooper- 
ative regions. Detecting amorphous excitations is not an 
easy task, as we lack a traditional order parameter (like 
the magnetization or the density in the magnetic or gas- 
liquid first order transitions) able to distinguish immedi- 
ately the presence of droplets of different phases. Hence, 
we need a protocol to artificially build amorphous exci- 
tations in the system. Note that the (low-temperature) 
excitations we are aiming to study here are different from 
the dynamically correlated regions (dynamical hetero- 
geneities) observed at higher temperatures, which are not 
necessarily associated with an energy cost related to their 
size and which may be string-shaped [1, 0| • 

The core of our idea to mimic the formation of amor- 
phous excitations is the following: we consider two inde- 
pendent equilibrium configurations a and (3 and simply 
exchange all the particles contained within a sphere of 
radius R. In this way we directly know the size and the 
position of the excitations and the excess energy cost dur- 
ing their evolution in time. The energy cost due to the 
formation of a droplet of phase a within a different phase 
P (and vice-versa) can then be studied. 

More in detail, this is what we do. We consider pairs 
of independently thermalized configurations, Ca and Cp, 
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FIG. 1: Construction of the mixed configurations Ca+p and 
C/s+a using configurations Ca and Cp. 



and from these we create a mixed configuration: all parti- 
cles within a sphere of fixed radius 7? of Ca are moved to a 
spherical cavity of the same shape and size in the config- 
uration Cf}] conversely, the particles within the sphere of 
Cf3 are moved in the spherical cavity of Ca as it is shown 
in figure [H In this way two new configurations arise, 
Ca+f3 and C/3+a- In each initial configuration the cavity 
is chosen in order to preserve the concentration of the 
two species of particles in the mixed configuration Ca+p- 
We can then study both the energetic and geometric evo- 
lution of these excitations at the same temperature as we 
used to obtain the initial configurations. 

By using this protocol, we know exactly two very im- 
portant things: 1. where the excitation is; 2. what its 
geometry is, and hence what coordinate is orthogonal to 
the interface between the two phases (in our case the 
radial one). Such knowledge is the only reason why we 
are able to compute and say something about these exci- 
tations. For spontaneous excitations, of course, we lack 
both pieces of information. 

This means, in particular, that we can measure the 
excess energy AEa/s of the excitation, due to the interface 
between a and 
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where Vij{ri — rj) is the pair potential, Eaf3 is the energy 
of the mixed configuration (originally a inside, /3 outside) 
and E™^ {E^^) is the energy of the particles inside (out- 
side) the sphere in configuration a (/3). As we shall see, 
all the relevant information about the exponents ■0 and 6 



comes from the excess energy, complemented by a study 
of the geometric properties of the excitation's surface. 

The system we consider is a binary mixture of soft par- 
ticles, a fragile glass-former i,^, i9|. In this model system, 
the particles are of unit mass and they belong to one of 
the two species 7 = 1,2, present in equal amount and 
interacting via a potential: 
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N 



V = ^y.,(|r,-r,|)=5] 



i<j 
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The radii are fixed by U2lo\ ~ 1.2 and setting the 
effective diameter to unity, that is (2(Ti)'^ -f 2(cri -I- a^Y + 
(2(72)'^ = 4Zq , where is the unit of length. The den- 
sity Ya p — N/V in units of , and we set Boltzmann's 
constant fc^ = 1. A long-range cut-off at Tc — ^/2> is 
imposed. The thermodynamic quantities of this system 
depend only on F = p/T^^'^, with T the temperature of 
the system [lO[. Here p ~ \ Iq^ thus the thermodynamic 
parameter F is F = 

The presence of two kinds of particles in the system 
strongly inhibits the crystallization and allows the ob- 
servation of the deeply supercooled phase. Moreover an 
efficient Monte Carlo (MC) algorithm (swap MC [TT]) is 
able to thermalize this system also at temperature below 
the Mode Coupling temperature (Fc = 1.45) [13] ■ 

In this work we study thermalized configurations of 
a system with N — 16384 particles confined in a pe- 
riodic box. The temperatures considered, correspond- 
ing to F = 1.49, 1.47, 1.44, 1.42, and 1.35, are two be- 
low and three above the Mode Coupling temperature T^. 
T ~ O.SgTc, O.gSTc, imr^, 1.09Tc, 1.33Tc. For each tem- 
perature we created mixed configurations as explained 
above from a collection of about 10 independent equi- 
librium configurations, using spheres of sizes between 3 
and 8 (in units of Iq). The equilibrium configurations 
were produced using the swap MC algorithm, but the 
relaxation of the mixed configurations was followed us- 
ing standard Metropolis MC. The upper limit of the size 
of the considered excitations is due to the emergence of 
boundary condition effects for larger droplets in the al- 
ready numerically challenging system with N — 16384 
particles. 



III. DETERMINATION OF tP AND 6 

As discussed in the introduction, the relaxation in 
deeply supercooled liquids proceeds through the acti- 
vated rearrangement of clusters of correlated particles, 
and to these rearrangments a barrier is associated, as- 
sumed to scale with a power of their size. As a result, 
the relaxation time is exponential in a power of the size of 
these regions, eq. ©. In this framework we expect that 
the timescales involved in the formation and in the relax- 
ation of the cooperative regions are in fact the same: each 
excitation relaxes through the cooperative rearrangement 
of new excitations. This means that we can follow the 
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process of relaxation of an artificially produced excita- 
tion, rather than detect the spontaneous formation of an 
excitation. 
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FIG. 2: {AE{R, t)) vs. t for R ranging from 4 to 8, left to right. 
Points are numerical data, curves are power law {AE{t, R)) ~ 
t ""'^^ fits. The fits are used to give an estimate of t(R), 
represented by crosses on the threshold line in figure. The 
values of 7 are within the range 0.2 < 7(-R) < 0.4. Inset; 
Full AE vs. t time series for R — 5, showing the kink that 
separates the fast decay and the slow relaxation regime. 
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FIG. 3: riog(r) vs. ii for T = O.SQTc; errors are evaluated 
using the bootstrap method. The linear behaviour on log-log 
plot allows for a scaling ansatz t{R) ~ exp{R/T). 

In practice, we study how the excess energy AE(t) 
of the excitation relaxes with time. The artificial con- 
struction of the excitations has a very large effect on the 
excess energy only for the first few MC steps (see Fig- 
ure [21 inset), so we disregard the data before the kink 
in the AE{t) curve. Figure [2] (obtained at T = 0.89Tc) 
clearly shows that the rate of relaxation of the excita- 
tions changes with their size R: bigger spheres relax over 
larger time scales. Simply by fixing a threshold value for 



AE{t), and measuring the time needed to drop below 
this value, we obtain an estimate of the time r needed to 
relax an excitation of size R. To actually find the time 
at which the threshold is reached, we need a way to in- 
terpolate between the data points at long times. This 
we do by fitting a power law to the data. We consider 
the range t > 100 MC steps. We perform the same pro- 
cedure at different sizes R of the sphere thus obtaining 
a function t{R). We have checked that the exponent in 
t{R) curve below is insensitive to changes in the thresh- 
old from AEth = 100 to AEth = 175. 

According to Tlog(r) has to scale as R^. In the 
log-log plot in figure [3] we report Tlog(r) vs. R for our 
lowest temperature. The data lie with good approxima- 
tion on a straight line, thus confirming that the process 
of relaxation of the excitations indeed follows the Arrhe- 
nius law. A fit of the exponent gives tp = 1.01 ± 0.04, 
and hence we conclude, 



iP = l. 



(14) 



When a region rearranges, it is natural to expect that 
some excess energy is stored at the interface between the 
new configuration and the old one. The exponent 9 reg- 
ulates how this interfacial energy scales with the size R 
of the excitation. 



AE = YR'^, 



(15) 



where the quantity Y is the (generalized) surface tension. 
As we have seen in section HI the exponent is crucial in 
order to discover the temperature dependence of the cor- 
relation length. Note however that 9 is an asymptotic 
exponent, and that subleading corrections to eq. 
are in general important for small sizes (see eq. be- 
low). These corrections are expected from curvature (as 
in liquid- liquid interfaces [l^l) or disorder effects (as in 
the random field Ising model [T^ . or the random bond 
Ising model [Tl|). 

Figure S] shows how the excess energy AE defined 
in p2)) scales with R for different temperatures and times. 
As we have seen, the excitations decay with time, and 
therefore for long times the dependence of AE on R be- 
comes rather hazy. However, for intermediate times and 
for all analyzed temperatures it is clear that a power law 
R^ seems to reproduce the asymptotic behavior rather 
satisfyingly. Attempting to fit the data at long times 
with a single power leads to 9 > 2, but 9 cannot grow 
with time. Indeed, the data at very short times follow 
an R^ law for all sizes (as it should by construction since 
the potential energy is short ranged). An exponent that 
grows with time is unacceptable, because it would im- 
ply that for sufhciently large i?, AE increases with time 
rather than decreasing. From this analysis we therefore 
conclude 

61:= 2. (16) 

The same value of 9 was found in [l^ using inherent 
structures. However, we remark that the present data are 
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FIG. 4: Time dependence of the excess energy AE{t) vs. R 
at four different temperatures. From left to right and top to 
bottom: T = O.SgTc, l.OSTc, 1.09Tc, LSSTc Time is measured 
in MonteCarlo steps. 



obtained at finite temperature, using equilibrium config- 
urations rather than potential energy minima. Hence, 
the value 9 = 2 seems to be quite robust. Moreover, the 
data also show that a nonzero surface tension Y survives 
for quite some time after the formation of the droplet, 
especially at the lowest temperatures. 



IV. CONSISTENCY WITH THE ROUGHENING 
EXPONENT 

Here we study the geometrical properties of the exci- 
tation interfaces, and find further support for = 2. 

In disordered systems interfaces are typically rough. 
The roughening of interfaces has been investi gate d at 
length since the directed polymer (DP) problem [15| and 
the Random Field or Random Bond Ising Model (RBIM) 
studies [TtI . [TH . The signature of roughening is the fact 
that the interface thickness w grows with the linear size 
of the interface itself. In our case it corresponds to the 
linear size R of the excitations: 



(17) 



where 7 is the so-called roughening exponent. We want 
to measure the value of 7 for the interfaces of the amor- 
phous excitations, and to link it to the exponent 6. How- 
ever, in order to study the roughening properties of the 
excitations we must use inherent structures (ISs), namely 
minima of the potential energy, rather than thermal con- 
figurations as wc have done up to now. The reason is 
that the roughening mechanism is ruled by a zero tem- 
perature fixed point, so that working at nonzero T would 
needlessly introduce the complication of treating thermal 
fluctuations. For the same reason, we focus on ISs ob- 
tained by equilibrium configurations at the lowest avail- 
able temperature, T = 0.89Tc. 



The procedure to create the excitations with the ISs is 
very similar to the one described in section [Hi We switch 
two spheres within two IS configurations C^^ and C^"^, 



to 



produce the mixed configuration Ca+/3- But such configu- 
ration is of course not an IS itself, so we must find the new 
minimum of the potential energy, C^^ . We do this by first 
performing 100 Monte Carlo steps at T = 0^ followed by 
an optimized quasi-Newton algorithm [limited memory 
Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) [H]. We 
used a collection of 16 ISs. For each pair of configura- 
tions we used spheres with radii between 1.5 and 8.5 in 
units of Iq. 




FIG. 5: Radial local overlap of a sphere of radius R = 8. 
For the estimate of w we look only at the outer part of the 
interface. This side is not affected by finite size effects in small 
spheres. 




^ This first part is required because a single pair of particles very 
close together produces a huge gradient that tends to destabilize 
the minimizer. 
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FIG. 6: Scaling of the interface thickness w with the size R 
of the sphere. 

In order to define the thickness w of the excitation 
we can look at the overlap q{r) (r is the distance from 
the origin) between the just-switched configuration Ca+f} 
and its relative minimum C^^ (for a definition of the local 
overlap see appendix]^ . The overlap is small close to the 
interface, where particles have moved most, while it gets 
close to one away from it. Hence, we can define w as the 
thickness of the region for which the local radial overlap 
of the excitation is smaller than an arbitrary threshold 
value qth- In figure [Sj the local overlap is plotted as a 
function of the distance r from the center of the sphere. 
Since in very small spheres the overlap has no room to 
reach high enough values in the inside, we actually define 
w as sketched in figure [5l In figure [6l we report w as 
a function of the radius R of the excitation in a log- 
log plot. We conclude that the excitations' interfaces 
roughen according to relation p?]) . with 

7 = 0.62 ±0.01. (18) 

A roughening exponent smaller than one implies that the 
ratio w / R between width and linear scale of the surface 
decreases for larger spheres. Thus large excitations have 
relatively thin interfaces. This is clearly shown in figure[71 
where two excitations with different radius (i? = 4 vs. 
i? = 8) are compared in a coordinate system where all 
lengths are rescaled by R, so that both rescaled spheres 
have virtual radius unity. The rescaling emphasizes the 
thick interface of the smaller excitation compared to the 
sharper interface of the larger excitation. 




FIG. 7: Local overlap qai3{r) between C^f} and Ca+fi for con- 
figurations with interface placed at 7? = 4 (left) and R = 8 
(right). Plots are rescaled so that both spheres appear as 
having the same size. Low overlap (dark grey) at interfaces 
indicate major rearrangements of particles, while far from in- 
terfaces the overlap with initial configurations is high (light 
grey). Clearly the larger sphere (right) has a relatively thin- 
ner, or smoother, interface. 

To understand why roughening occurs we can think 
about the Random Bond Ising Model (RBIM). This is 



an Ising spin model where the nearest neighbours bonds 
are random, albeit typically positive. In such a system 
the position of a domain wall strongly depends on the 
disorder, since the weak bonds are more likely to be bro- 
ken. On one hand, a smooth domain wall is preferable, 
as it would break the smallest number of bonds. On 
the other hand, some suitable deviation from smoothness 
could induce the breaking of weaker bonds and hence a 
lower energy cost. Hence, a rough interface is the re- 
sult of a complicated optimization problem: the cost of 
a large number of broken bonds is balanced by the gain 
due to the presence of very weak bonds among them. As 
a result, in a disordered system a rough interface can be 
energetically favoured with respect to a smooth interface. 

The interesting point is that in the context of elastic 
manifolds in random media ^20l . [2]| , a precise relation ex- 
ists between the energy gain we just mentioned and the 
roughening exponent 7, and in such relation the expo- 
nent 6 comes into play. The energy gain by roughening 
appears as a negative correction to the ground state en- 
ergy of the manifold, 

Es^YR^ - BR^'^ , d = 3. (19) 

Our hypothesis is that the interface energy AEafs of the 
amorphous excitations can be described as in (|19p by a 
leading term (due to a generalized surface tension) plus 
a sub-leading correction due to roughening. In random 
manifolds 9 = d—1, whereas we do not know 9. We thus 
try to find 9 by taking the value of 7 found in (fT5)) and 
using eq. (|19p to fit AEap with as a fitting parameter. 
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FIG. 8: (AEi^R)) behaviour with 7? at different temper- 
atures T ranging from O.SgTc to L33Tc. The log-log plot 
enhances the power law trend at large R. Lines are fits to 
eq. (|19p with and 7 fixed (see text). 

The interface energy AEap is calculated as in eqs. 
except that Eafj is the energy of C^^. In figure [8] we plot 
AEap vs. R at different temperatures. A deviation from 
a simple power law at small sizes is evident. A fit of 
AEajs at the lowest temperature T = 0.89Tc with the 
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functional form ((T9)) (not shown) and 7 = 0.62 gives 

61 = 2.08 ± 0.04. (20) 

Reversing the procedure (i.e. fixing 9 = 2 and fitting 7 in 
using (HH])) gives 7 = 0.75. Finally, fixing both 9 = 2 and 
7 = 0.62 and 9 = 2, still gives an excellent fit of AEap vs. 
i? (figure [8]). Since we know that 9 cannot exceed 2, this 
consistency check confirms the value of 9 found before. 

V. DISCUSSION 

Our result 9 = 2 is somewhat sensible and not par- 
ticularly exciting: it is basically telling us that disorder 
in a supercooled liquid is not strong enough to change 
in any exotic way the leading term of the surface energy 
cost: surfaces remain surfaces, albeit a bit rough. 9 is 
not actually a part of AG theory, so it is more proper to 
compare with RFOT, where 9 was originally introduced. 
It must be said that the value 0rfot — d/2 derived in 
[31 using RG arguments always had its greatest appeal 
in the fact that, together with -iARFOT = ^rfot, it gave 
back the VFT equation (fTTjl . In fact, the arguments used 
in Q to fix do not belong to RFOT itself, and other 
values are in principle compatible with the conceptual 
structure of RFOT. On the contrary, the value found for 
ijj within RFOT has to do with a crucial aspect of the the- 
ory. Hence, it seems that the real interesting comparison 
is about the exponent ip, after all. 

The value ^p = 1 we find implies that the barrier for 
the rearrangement of a correlated region scales linearly 
with its size, 

A^e- (21) 

The proposal of AG, A ~ has always seemed a bit 
exaggerated, since one can imagine several ways for the 
system to pay less than the entire volume to rearrange a 
region. One of this ways, and not the most pedestrian, is 
suggested by the fact that the greatest part of the energy 
necessary to create the excitation is stored in the inter- 
face. Hence, one may expect that the barrier scales with 
the same exponent as the surface energy cost, i.e. 9. This 
is, basically, the idea of RFOT. We note that such idea is 
deeply rooted in the theory of nucleation, which was a in- 
deed a source of inspiration for the original formulation 
of RFOT. A proportional relationship between inverse 
diffusion constant and the exponential of the number of 
particles belonging to a dynamically correlated cluster 
has been reported for a model of water [22], but this is 
not necessarily in contradiction with the above. The dy- 
namical clusters and the excitations of RFOT (or AG's 
cooperatively rearranging regions) are probably different 
entities, as mentioned above. Furthermore, the simu- 
lations of ref. Il^ correspond to a different temperature 
regime (above the mode-coupling temperature) and the 
objects found in that work are rather noncompact, with 
different volume/surface ratio from the excitations we are 
studying. 



According to RFOT there are two competing forces: 
the free-energy cost to create an excitation, scaling as 
YR^, and the configurational entropy gain due to the 
change of state of the rearranging region, scaling as 
TScR''^. In these two expressions Y is the surface ten- 
sion and Sc is the configurational entropy. Hence, the 
total free energy for the formation of the excitation is, 
according to RFOT, 

AF{R) = YR'^ - TScR'^. (22) 

At this point, RFOT, in perfect analogy with nucleation 
theory, proceeds by finding the maximum of such non- 
monotonous function (recall that 9 < d). This maximum 
provides two essential pieces of information: First, the 
position of the maximum, R = ^, gives the critical size of 
the rearranging region, i.e. the mosaic correlation length 
(cf. eq.[3D, 

Second, and most important for us now, the height of the 
maximum, AF{R = ^), gives the size of the free-energy 
barrier to be crossed to rearrange the region, 

A = AF{R = C)- (24) 

This fixes V'RFOT — ^rfot, and it coincides with the 
intuitive notion that the barrier should scale the same as 
the interface cost. 

This last result, however, is at variance with what we 
find here, relation (|2ip : in fact, whatever one thinks 
about our numerical result for 9, a value of 9 as small 
as 1 seems rather unlikely. In any case, it is important to 
emphasize that "tjj = 9 is a, consequence of the maximiza- 
tion of eq. (|22p . which in turn follows from the nucleation 
paradigm. It has been noted, however, that nucleation 
is perhaps not a fully correct paradigm to describe the 
formation of amorphous excitations within a deeply su- 
percooled liquid [23|. The essence of RFOT, namely the 
competition between a surface energetic term and a bulk 
entropic term, retains its deepest value even if we do not 
cast it within the strict boundaries of nucleation theory. 
The value of the correlation length may come from the 
point where the two contributions balance, rather than 
from the maximum of , and (for obvious dimensional 
reasons) one gets the same expression (^5)) for ^ (up to 
an irrelevant constant), while ip remains undetermined. 
These points are discussed in depth in 23]. Here we 
simply note that our present results are quite compatible 
with RFOT in the form it has been recast in [2^ , without 
reference to a nucleation mechanism. 

Regarding the comparison between 9 and "0 there is a 
final point we have to discuss. The reader familiar with 
spin-glass physics will probably remember the Fisher- 
Huse (FH) inequality, |24| . 

^>9, (25) 
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which is plainly violated by the values we find here. The 
physical motivation of (PS)) is basically the following: if 
we represent the excitation as an asymmetric one dimen- 
sional double well, where the abscissa is the order param- 
eter and the ordinate is the energy of the excitation, the 
height of the barrier (which scales as £,^) is always larger 
than (or equal to) the height of the secondary minimum 
(which scales as ^^), and hence ^ > 0. 

How comes, then, that we find ip < 91 The FH argu- 
ment was formulated in the context of the droplet pic- 
ture for spin-glasses, where there are only two possible 
ground states. In supercooled liquids, in contrast, a re- 
arranging region can choose among an exponentially large 
number of target configurations. Under this conditions, 
even though the FH bound still applies to the energy 
barrier, there may be a nontrivial entropic contribution 
that decreases the free energy barrier to rearrangement. 
We determine if) by measuring a time, and hence a free- 
energy barrier, not an energy barrier, so that the FH 
constraint does not necessarily hold in the case of su- 
percooled liquids. This entropic effect is absent in the 
original FH argument due to the lack of exponential de- 
generacy of the target configurations (in fact, one would 
expect relation ([25]) to hold even in the mean-field pic- 
ture of spin-glasses, where the number of ground states is 
large but the configurational entropy is still zero). How 
the FH argument should be modified in the presence of 
such large entropic contribution is however not clear at 
this point. 

VI. CONCLUSIONS 

We have determined numerically (in d — i) the expo- 
nents linking the size of rearranging regions to barrier 
height (V') and to surface energy cost [9). This is to our 
knowledge the first direct measure of these exponents. 
We find 

ip = l, 9 = 2. (26) 

Both values are in disagreement with those assumed by 
the AG and RFOT schemes. However, if we stick to 
eq. ([5]), we still obtain for the relaxation time the VFT 
relation pTjl . It seems that, albeit changing all cards 
on the table, we managed to get back the most used fit- 
ting relation in the physics of glass- forming systems. We 
stress that there is no particular reason to stick to such 
VFT In fact, as it has been remarked many times before, 
a generalised VFT with an extra fitting exponent v, such 
that logr ^ (T — Tk)~'^ , would do an even better job in 
fitting the data. Yet, to get back VFT as the product of 
the independent numerical determination of two rather 
different exponents, remains a rewarding result to some 
extent. 



Finally, let us remark that the disagreement between 
our exponents and those proposed in the original RFOT 
do not imply as harsh a blow to RFOT as it might seem 
at first. RFOT can be cast in a form [23j that retains its 
most essential aspect, namely the competition between 
a surface energy cost and a bulk energy gain, without 
using nucleation theory. If one does this, the exponents ip 
and 9 remain unrelated and compatible with our findings. 
Within this context, our result ip = 1 seems to be an 
indication that RFOT is a better theory if one does not 
push the nucleation analogy too hard. 
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APPENDIX A: DEFINITION OF THE LOCAL 
OVERLAP 

A suitable definition of the overlap, for the off-lattice 
system considered, is given using a method similar to 
that used in [2^ : we divide the system in 64'^ small cubic 
boxes with side L/64 and, having two configurations a 
and r, we compute the quantity 

Qarix, y, z) = n„{x, y, z)nr{x, y, z) , (Al) 

where n„{x,y, z) is 1 when the box with center at coor- 
dinates x,y,z contains at least one particle and it is 
when the same box is empty. To each box we assign the 
weight 

n^{x,y,z) +nr{x,y,z) 
Wcr{x,y,z) = ^ . (A2) 

The global overlap q^r between these two configurations 
a and r is given by 

^crri^i: yij Zi)'WaT{Xi: yii Zi) f \ o\ 

laT — J- N • y^'J) 

where i runs over all boxes in the system. 

We are interested in the local value of the overlap 
9(1/3 (^) distance r from the centre of the sphere. The 
definition of the local overlap qafsir) for a spherical 
corona between r — dr and r + dr is given by consid- 
ering in (|A3p only the sum of boxes belonging to this 
region. 
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